Existence and uniqueness for a nonlinear 
parabolic/Hamilton- Jacobi coupled system describing 
the dynamics of dislocation densities 



Hassan IBRAHIM* 
February 2, 2008 



Abstract 

We study a mathematical model describing the dynamics of dislocation densities in crystals. This 
model is expressed as a one-dimensional system of a parabolic equation and a first order Hamilton- 
Jacobi equation that are coupled together. We show the existence and uniqueness of a viscosity solution 
among those assuming a lower-bound on their gradient for all time including the initial data. Moreover, 
we show the existence of a viscosity solution when we have no such restriction on the initial data. We 
also state a result of existence and uniqueness of an entropy solution of the system obtained by spatial 
derivation. The uniqueness of this entropy solution holds in the class of "bounded from below" solutions. 
In order to prove these results, we use a relation between scalar conservation laws and Hamilton- Jacobi 
equations, mainly to get some gradient estimates. This study will take place in M, and on a bounded 
domain with suitable boundary conditions. 

Resume 

Nous etudions un modele mathematique decrivant la dynamique de densites de dislocations dans les 
cristaux. Ce modele s'ecrit comme un systeme ID couplant une equation parabolique et une equation 
de Hamilton- Jacobi du premier ordre. On montre I'existence et I'unicite d'une solution de viscosite dans 
la classe des fonctions ayant un gradient minore pour tout temps ainsi qu'au temps initial. De plus, on 
montre I'existence d'une solution de viscosite sans cette condition sur la donnee initiale. On presente 
egalement un resultat d'existence et d'unicite pour une solution entropique d'un systeme obtenu par 
derivation spatiale. L'unicite de cette solution entropique a lieu dans la classe des solutions minorees. 
Pour montrer ces resultats, on utilise une relation entre les lois de conservation scalaire et les equations 
de Hamilton- Jacobi, principalement pour obtenir des controles du gradient. Cette etude a lieu dans M 
et dans un domaine borne avec des conditions aux bords appropriees. 
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1 Introduction 



1.1 Physical motivation 

A dislocation is a defect, or irregularity within a crystal structure that can be observed by 
electron microscopy. The theory was originally developed by Vito Volterra in 1905. Disloca- 
tions are a non-stationary phenomena and their motion is the main explanation of the plastic 
deformation in metallic crystals (see |28ill9j for a recent and mathematical presentation). 

Geometrically, each dislocation is characterized by a physical quantity called the Burgers 
vector, which is responsible for its orientation and magnitude. Dislocations are classified as 
being positive or negative due to the orientation of its Burgers vector, and they can move in 
certain crystallographic directions. 

Starting from the motion of individual dislocations, a continuum description can be derived 
by adopting a formulation of dislocation dynamics in terms of appropriately defined disloca- 
tion densities, namely the density of positive and negative dislocations. In this paper we are 
interested in the model described by Groma, Csikor and Zaiser [18], that sheads light on the 
evolution of the dynamics of the "two type" densities of a system of straight parallel dislocations, 
taking into consideration the infiuence of the short range dislocation-dislocation interactions. 
The model was originally presented in X (0, T) as follows: 



de+ , d 



r 



[Tsc + Teff) - AD- 



b d 



— -h — 

dt dr 



{e+ + e-) dr 



0, 

(1.1) 

0. 



Where T > 0, r = (x,y) represents the spatial variable, b is the burger's vector, 9~^{Y,t) 
and 6~{r,t) denote the densities of the positive and negative dislocations respectively. The 
quantity A is defined by the formula A = ///[27r(l — z^)], where /i is the shear modulus and u 
is the Poisson ratio. -D is a non-dimensional constant. Stress fields are represented through 
the self-consistent stress Tsc{r,t), and the efi'ective stress re//(r,i). ^ denotes the gradient 
with respect to the coordinate vector r. An earlier investigation of the continuum description 
of the dynamics of dislocation densities has been done in [T^. However, a major drawback 
of these investigations is that the short range dislocation-dislocation correlations have been 
neglected and dislocation-dislocation interactions were described only by the long-range term 
which is the self-consistent stress field. Moreover, for the model described in [17], we refer the 
reader to [HI [12] for a one-dimensional mathematical and numerical study, and to [4] for a 
two-dimensional existence result. 

In our work, we are interested in a particular setting of (jl.ip where we make the following 
assumptions: 

(al) the quantities in equations p.ljl are independent of y, 
(a2) b = (1, 0), and the constants A and D are set to be 1, 
(a3) the effective stress is assumed to be zero. 

Remark 1.1 (al) gives that the self-consistent stress Tgc is null; this is a consequence of the 
definition of Tsc (see fl^)- 
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Assumptions (al)-(a2)-(a3) permit rewriting the original model as a ID problem in M x (0,T): 



'ix,t)- (0+ix,t) 



6t{x,t)-9-{x,t) 

e+{x,t) + e-{x,t) 

« / N f^, ^ ( Qt{x,t) - 0~(x,t) 

We consider an integrated form of (|1.2p and we let: 



0, 



(1.2) 



p^ = 0^^ = 9+ + e', p = p+-p- and K = p+ + p", 



(1.3) 



in order to obtain, for special values of the constants of integration, the following system of 
PDEs in terms of p and k, : 



i^tUx = PtPx in Qt 

.0/ 



X (o,r), 



k(x,0) = K^(x) 



m 



and 



Pt = Pxx in Qt, 

p(x,0) =p°(x) 



m 



(1.4) 



(1.5) 



where T > is a fixed constant. Enough regularity on the initial data will be given in order to 
impose the physically relevant condition, 



> \pI\ . (1.6) 

This condition is natural: it indicates nothing but the positivity of the dislocation densities 
0^(x,O) at the initial time (see (II. 3p ). 



1.2 Main results 

In this paper, we show the existence and uniqueness of a viscosity solution n of (II. 4p in the 
class of all Lipschitz continuousviscosity solutions having special "bounded from below" spatial 
gradients. However, we show the existence of a Lipschitz continuousviscosity solution of (|1.4p 
when this restriction is relaxed. A relation between scalar conservation laws and Hamilton- 
Jacobi equations will be exploited to get almost all our gradient controls of n. This relation, 
that will be made precise later, will also lead to a result of existence and uniqueness of a 
bounded entropy solution of the following equation: 

{a ( PxPxx \ . ^ 

''=\-^)x (1.7) 
e{x,^)=d^{x) in M, 

which is deduced formally by taking a spatial derivation of p.4p . The uniqueness of this 
entropy solution is always restricted to the class of bounded entropy solutions with a special 
lower-bound. 

Let Lip{W) denotes: 

Lip(M) = {/:Mh^M; /is a Lipschitz continuous function}. 

We prove the following theorems: 
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Theorem 1.2 (Existence and uniqueness of a viscosity solution) 

Let T > 0. Take G Lip{R) and G C^(]R) as initial data that satisfy: 

K° > ViPx)'^ + a.e. m M, (1.8) 

for some constant e > 0. Then, given the solution p of ^-5]) . there exists a viscosity solution 
K G Lip{QT) of ( f j.^D , unique among the viscosity solutions satisfying: 

Theorem 1.3 (Existence and uniqueness of an entropy solution) 

LetT>0. Take 6*° G L°°(R) and p° G C^(M) smc/i that, 

0^ > ViPx)^ + m M, 

/or some constant e > 0. T/ien, ^/lere exists an entropy solution 9 G L^{Qt) of 111-7]) . unique 
among the entropy solutions satisfying: 

6 > \/ Px + in Qt- 

Moreover, we have 9 = Kx, where k is the solution given by Theorem 

The notion of viscosity solutions and entropy solutions will be recalled in Section [2l We now 
relate these results to our one-dimensional problem (|1.2p . Remarking that Px = 9~^ — 9~ and 

Kx = 9^ + 9~ , we have as a consequence: 

Corollary 1.4 (Existence and uniqueness for problem (II. 2j) ) 

Let T > 0. Let 9^ and 9q he two given functions representing the initial positive and negative 
dislocation densities respectively. If the following conditions are satisfied: 

(1) 9+ - 9, G Co-(M), 

(2) 9+, 9o G L-(M), 
together with, 

Ot + 9o> \J {9t - OqY + e"^ a.e. in M, 

then there exists a solution {9~^,9~) G {L°°{Qt))'^ to the system in the sense of Theorems 

Ij.^l and \l.^ unique among those satisfying: 

9+ + 9~ > y/{9+ - 0-)2 + e2 a.e. in Qt- 

Remark 1.5 Conditions (1) and (2) are sufficient requirements for the compatibility with the 
regularity of p^ and previously stated. 

Theorem 1.6 (Existence of a viscosity solution, case e = 0) 

Let T > 0, K° G Lip{R) and p° G C^(M). // the condition / flTl) is satisfied a.e. in M, then 
there exists a viscosity solution k G Lip{Q'r) of |.?.^[ ) satisfying: 

i^x ^ \Px\ o,.e. in Qt- (1-9) 
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Remark 1.7 In the limit case where e = 0, we remark that having U . 9\} was intuitively expected 
due to the positivity of the dislocation densities 9^ and 9~ . This reflects in some way the well- 
posedness of the model of the dynamics of dislocation densities. We also remark that our 

result of existence of a solution of (T^ under lll.9\) still holds if we start with = = on 
some interval of the real line. In other words, we can imagine that we start with the probability 
of the formation of no dislocation zones. 

Problem with boundary conditions. 

We consider once again problem (|1.4I) . similar results to that announced above will be shown 
on a bounded interval of the real line with Dirichlet boundary conditions (see Section (S]). This 
problem corresponds physically to the study of the dynamics of dislocation densities in a part 
of a material with the geometry of a slab (see |18|). 

1.3 Organization of the paper 

The paper is organized as follows. In Section [2l we start by stating the definition of viscosity 
and entropy solutions with some of their properties. In Section [3l we prove the existence and 
uniqueness of a viscosity solution to an approximated problem of (|1.4p . namely Proposition 
13.11 and we move on, giving additional properties of our approximated solution (Proposition 
13. 2p and consequently proving Theorems 11.21 and 11.31 In Section HJ we present the proof of 
Theorem 11.61 Section [5] is devoted to the study of problem (|1.4I) on a bounded domain with 
suitable boundary conditions. Finally, Section[6]is an appendix containing a sketch of the proof 
to the classical comparison principle of scalar conservation laws adapted to our equation with 
low regularity. 

2 Notations and Preliminaries 

We first fix some notations. If O is an open subset of M", k is a positive integer, we denote 
by C'^(O) the space of all real valued k times continuously differentiable functions. Cq{^) is 
the subspace of C^{Q) consisting of function of compact support in and C^{Q) = C'^{Q) H 
W'^'°°{i}) where VF^'°°(0) is defined below. Furthermore, let UC{0,) and Lip{^l) denote the 
spaces of uniformly continuous functions and Lipschitz continuousfunctions on respectively. 
The sobolev space W^'P{VI) with m > 1 an integer and p : 1 < p < oo a real, is defined by 

Va with |a| < n 3fa G L^(r2) such that ^ 

/ uD^cj) = (-1)1"! / V0 G Co^(O) ( ' 
Jn Jn ^ 

where we denote D"u = fa- This space equipped with the norm 

0<|o|<n 

is a Banach space. In what follows, T > 0. A map m : [0, oo) i— > [0, oo) that satisfy 

• m is continuous and non-decreasing; 

• lim m(x) = 0; 
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• m{a + b) < m{a) + m{h) for a, 6 > 0; 

is said to be "a modulus", and UCxiyt. x [0, T]) denotes the space of those u G C(r2 x [0, T]) for 
which there is a modulus m and r > such that 

\u{x^ t) — u{y, t)\ < m{\x — y\) for x,y £ ^l, \x — y\ < r and t £ [0, T]. 

We will deal with two types of equations: 

1. Hamilton- Jacobi equation: 

Ut + F{x,t,Ux) 
u{x, 0) = u^{x) 

2. Scalar conservation laws: 

Vt + {F{x,t,v))x 
v{x, 0) = v^{x) 

where 

F : M X [0, T] X M 
{x,t, u) 

is called the Hamiltonian in the Hamilton-Jacobi equations and the flux function in the scalar 
conservation laws. We will agree on the continuity of this function, while additional and specific 
regularity will be given when it is needed. 

Remark 2.1 We will use the function F as a notation for the Hamiltonian /flux function. 
Although F might differ from one equation to another, it will be clarified in all what follows. 

Remark 2.2 The major part of this work concerns a Hamiltonian/ flux function of a special 
form, namely: 

F{x,t,u) = g{x,t)f{u), (2.3) 

where such forms often arise in problems of physical interest including traffic flow ]31^ and 
two-phase flow in porous media /7^ . 

We start by defining the notion of viscosity solution to Hamilton-Jacobi equations (12.11) . 
and entropy solution to scalar conservation laws (|2.2p with a flux function given by Remark [221 
as well as some results about existence, uniqueness, and regularity properties of these solutions. 
We will end by a classical relation between these two problems. These results will be needed 
throughout this paper, precise references for the proofs will be mentioned later on. 

2.1 Viscosity solution: definition and properties 

Definition 2.3 (|10|. Viscosity solution: non-stationary case) 

1) A function u G C{Qt]^) is a viscosity sub-solution of 

ut + F{x,t,Ux) = in Qt, (2.4) 



in Qt, 

in M, 

= in Qt, 

in M, 



F{x, t, u) 



(2.1) 



(2.2) 
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if for every 4> G C^{Qt), whenever u — (j) attains a local maximum at (xo,to) £ Qt, then 

4>t{xo,to) + F{xo,to,(pxixo,to)) < 0. 

2) A function u G C{Qt','^) is a viscosity super- solution of {2.4^ if for every (j) G C^{Qt), 
whenever u — <f> attains a local minimum at {xo,to) G Qt, then 

4>t{xo,to) + F{xo,to,(pxixo,to)) > 0. 

3) A function u G C{Qt',^) is a viscosity solution of {2. 4^ if it is both a viscosity sub- and 
super-solution of \2.4^ . 

4) A function u G C{Qt',^ is a viscosity solution of the initial value problem (2. 1]) if u is a 
viscosity solution of (2.4\ ) and u{x, 0) = n°(x) m K. 

It is worth mentioning here that if a viscosity solution of a Hamilton-Jacobi equation is dif- 
ferentiable at a certain point, then it solves the equation there (see [TOl Corollary 1.6]). An 
equivalent definition depending on the sub- and super-differential of a continuous function is 
now presented. This definition will be used for the demonstration of Proposition 12.101 Let us 
recall that the sub- and the super-differential of a continuous function u G C(M" x (0, T)), at 
a point {x,t) G x (0, T), are defined as the closed convex sets: 



D^-uix,t) = G M" X 

liminf 

{y,s)^{x,t) \y — x\ -\- \s — t 



ujy, s) - u{x, t) - {p - {y - x) + a - {s - t)) ^ ^| 



and 



D^'+u{x,t) = {(p, 



a) G 



lij^g^p g) - ujx, t)-{p-{y-x)+a-{s- t)) ^ ^] 

{y,s)^{x,t) \y-x\ + \s-t\ ~ J 



respectively. 



Definition 2.4 (Equivalent definition of viscosity solution) 

1) A function u G C(IR" x (0,T)) is a viscosity super- solution of ^2. 1\} if and only if, for every 
{x,t) G M" X (0,r); 

y{p,a) £ D^'-u{x,t), a + F{x,t,p) >0. (2.5) 

2) A function u G C(M" x (0,T)) is a viscosity sub-solution of 12. 1\} if and only if, for every 
{x,t) G X (0,r); 

y{p,a) £ D^'+u{x,t), a + F{x,t,p) <0. (2.6) 

This definition is more local, for it permits verification that a given explicit function is a viscosity 
solution in a more classical way, i.e. using the derivative calculus. A similar definition, that 
will be used later, could be given in the stationary case. Let Q C be an open domain, and 
consider the PDE 

F{x,u{x),Vu{x)) =0, yxen, (2.7) 
where F : x M x M" M is a continuous mapping. 
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Definition 2.5 (Viscosity solution: stationary case) 

A continuous function u : M is a viscosity sub-solution of the PDE \2. 7\) if for any 

continuously differentiable function (j) : ^ M. and any local maximum xq S f2 of u — (j), one 
has 

F{xo,u{xo),V^{xo)) < 0. 
Similarly, if at any local minimum point xq £ il. of u — (j), one has 

F{xo,u{xq),V4){xq)) > 0, 

then u is a viscosity super- solution. Finally, if u is both a viscosity sub-solution and a viscosity 
super-solution, then u is called a viscosity solution. 

In fact, this definition is used for interpreting solutions of (|1.4p in the viscosity sense. Further- 
more, we say that u is a viscosity solution of the Dirichlet problem (|2.7p with u = C G C{dQ) 
if: 

(1) ueC{n), 

(2) u is a viscosity solution of (12. 7|) in fi, 

(3) n = C on dQ. 

For a better understanding of the viscosity interpretation of boundary conditions of Hamilton- 
Jacobi equations, we refer the reader to [21 Section 4.2]. 

Now, we will proceed by giving the main results concerning viscosity solutions of (|2.1I1 . In 
order to have existence and uniqueness, the Hamiltonian F will be restricted by the following 
conditions : 

(FO) F G C{R X [0,r] X M); 

(Fl) for each R > there is a constant Cr such that for all {x,t,p), {y,t,q) G R x [0,r] x 
[-R,R], 

I F{x,t,p) - F{y,t,q) \ < Cr{ \p-q\ + \x - y\)] 
(F2) there is a constant Cp such that for all G [0,T] x M and all a;, y G M, 

I F{x,t,p) - F{y, t,p)\< Cf\x -y\{l + \p\). 

We use these conditions to write down some results on viscosity solutions. 

Theorem 2.6 (Comparison, [9l Theorem 1]) 

Let F satisfy (F0)-(F1)-(F2). Ifu, n G UCx{Qt) o,re two viscosity sub- and super-solution of 
the Hamilton-, Jacobi equation \2. 1\) respectively, with 

u(x, 0) < ^(a;, 0) in M, 

then u < u in Qt- 

Theorem 2.7 (Existence, [9', Theorem 1]) 

Let F satisfy {F0)-{F1)-{F2). Ifu^ G UC(R), then |lJp has a viscosity solution u e UC.j:{Qt). 

Remark 2.8 The "comparison" theorem stated above gives the uniqueness of the viscosity so- 
lution. 
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Remark 2.9 In the case where the Hamiltonian has the form 



F{x,t,u) = g{x,t)f{u) 



the following conditions: 



(VO) / G Cl 



(M: 



R) 



(V2) g., G L-(Qt) 



imply (F0)-(F1)-(F2) together with the boundedness of the Hamiltonian. 

The next proposition reflects the behavior of viscosity solutions under additional regularity 
assumptions on and F. 

Proposition 2.10 (Additional regularity of the viscosity solution) 

Let F = gf satisfy (V0)-(V1)-(V2) . If G Lip{R) and u G UCx{Qt) is the unique viscosity 
solution of (2. 1]) . then u G Lip{QT). 

Proof. Consider the function defined on M x [0, T] by: 



By [20l Theorem 3], the function u satisfies, 

\u{x,t)\ < c*{\x\ + 1) for {x,t) G M X [0,r], 

where c and c* are two positive constants. Therefore, is a sublinear function for every 
time t G [0,T]. The function is defined via a supremum which is attained because of the 
sublinearity of the function u (a quadratic function always control a linear one); the supremum 
can be achieved at several points; let x^ be one of them, so we can write 




u^{x, t) 



u{xs,t) - e' 




We are going to prove that for {p, a) G M x M, we have: 




(2.8) 



Since {p,a) G D^'^u^{x,t), then we can write for {y,s) ~ {x,t) that. 



L = u'^{y, s) < u^{x, t) + a{s — t)+ p{y — x) + o{\s — t\ + \y — x\) = R, 



(2.9) 



where the left side L of (|2.9p satisfies 



L > u[z, s) — e' 



ks 



z-y 



z G M, 



(2.10) 



2e 
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and the right side R of (12. 9|) satisfies, 

I |2 

R < uixe,t) - e^*^^^ + a{s - t) + p{y - x) + o{\s -t\ + \y- x\). (2.11) 
Choose z such that z — y = x^ — x, then 

z = Xe + (y — x) ~ Xg, since y ~ x. (2-12) 
Combining (l2J|), ([2l^ . ([2llD and (l2T2]l together, we get 

/ / \ \ ks 1"^ I 

u(x£ + (y - x),s) - e < 



■u(x£, t) - e^* '^ + a(s - t) + ^(z - Xe) + o(|s - i| + |z - x^l). 



|2 



and hence, 

^(z, s) < u{xe,t) + (e"'' - e" 



2e 

+ a(s-t)+p(z-X£) + o(|s-t| + |2;-X£|). (2.13) 

We have 



(g.. _ ^ fce^* '"" J^' (g -t) + o{\s - t\), 



then using inequahty (|2.13p . we get 



u{z, s) < u{xs,t) + + fce^* J (s - t) 

+ ji(z - Xg) + o(|s - i| + |z - Xel), 

which proves that 

(a + fce^^^^^^.p) Gl)i'+n(x„t), 
and hence statement (|2.8I1 is true. Since u is a viscosity sub-solution of (|2.1I) . we have 



a + fce^* '^ Z'^' +F(xs,t,j)) < 0. 



2 



We use condition (Fl) with p = g, to get 



a + fee^*^^^ + F(x,t,p) < F{x,t,p)-F{xe,t,p), 

< Clx-Xel, 



2 



therefore, 



a + F(x,t,p) < C\x — Xs\ — ke 

< Cr, - k^, 

2e 

( kr^ 

< sup or 



kt\^ ^el 



2e 



r>0 
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2e / ' 



where = |x — x^l. At the maximum f, we have C = By choosing A; = we get 

a + F{x, t,p) < e. 

This inequahty shows that = — et is a viscosity sub-solution of (12.11) with v^{x,0) 
u^{x,0). By the comparison principle, we have 

v'^{x,t) — u{x,t) < sup(f^(a;, 0) — n°(x)), 

< sup{u^x,0) -u°{x)), 

|2 ■ 



< sup ( sup|n°(y) - ^1 ) , 
xeR \yeR I ^£ J / 



(I X — y\'^ 
7|x — y\ 

where 7 is the Lipschitz constant of the function u^, and r = \x — y\. This altogether shows 
the following inequality for x, y S M: 

I |2 2 

t) - e^*^^-^ < n^(x, t) < u{x, t)+et + ^. (2.14) 
Remark here that fc is a fixed; previously chosen constant. Inequality (|2.14p yields: 

I |2 / 2 \ 

u{y, t) - u{x, t) < e^*^^ + (^t + ^je = C/e + f3e, (2.15) 

where C = e^* and /3 = + . We minimize inequality (|2.15p over e to obtain, 

uiy,t)-u{x,t) < 2v^, 

< e^V2^t + ^\x-y\. 

Since this inequality holds Vx, y G M, exchanging x with y yields, 

\u{x,t) - u{y,t)\ < C{F,uo)\x - y\ Vx,yGM and iG[0,r]. 

This shows that the function u is Lipschitz continuous in x, uniformly in time t. To prove the 
Lipschitz continuity in time, we mainly use the result of [20j Theorem 3]) with the fact that 
ut = —F{x,t,Ux)^ and the boundedness of the Hamiltonian. □ 

Remark 2.11 It is worth mentioning that the space Lipschitz constant of the function u de- 
pends on C , where C appears in (Fl) forp = q, and on the Lipschitz constant 7 of the function 
uq. While the time Lipschitz constant depends on the hound of the Hamiltonian. 
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2.2 Entropy solution: definition and properties 
Definition 2.12 (Entropy sub-/super-solution) 

Let F{x,t,v) = g{x,t)f{v) with g, G L'^^{Qt;'M.) and f G C^(M;M). A function v £ 
L°°(Qt;IK) is an entropy sub-solution of (KM '^^^^ hounded initial data G L°°(M) if it 
satisfies: 

'ni{v{x,t))4>t{x,t) + <^{v{x,t))g{x,t)4>:r{x,t)+ 

f ^^'^^^ 

h{v{x,t))gx{x,t)(f){x,t) dxdt + I 'qi{v^{x))(f){x,0)dx > 0, 



V(/) G Co(M X [0,r);M+), for any non- decreasing convex function r]i G C^(M;M), $ G C^(R;M) 
such that: 

^' = f'vi, and h = ^-fri,. (2.17) 

An entropy super- solution of S2.S\} is defined by replacing in S2.16\) rji with r]d; a non-increasing 
convex function. An entropy solution is defined as being both entropy sub- and super-solution. 
In other words, it verifies 112. 16\) for any convex function r] G C^(M;M). 

A well know characterization of the entropy solution is that: 

Proposition 2.13 A function v G L°°{Qt) is an entropy sub-solution of if and only if 

V/c G M, G C^{R X [0,r);R+), one has: 

{v{x,t) - k)+(l)t{x,t) + sgn+{v{x,t) - k){f{v{x,t)) - f{k))g{x,t)Mx,t)- 



sgn^{v{x,t) -k)f{k)gx{x,t)<j){x,t) dxdt + / {v^ {x) - k)^ <j){x,Q)dx >Q, (2.18 



IS 



0/ 

u 

R 

Where = \{\a\ ± a) and sgn^{x) = ^{sgn{x) ± 1). An entropy super- solution of (KM) 
defined replacing in (2.18\} {■)~^ , sgn'^ by {•)~ , sgn~ . 

This characterization can be deduced from (|2.16|) . by using regularizations of the function 
(• — k)~^ . Also (|2.16|) may be obtained from (|2.18|) by approximating any non-decreasing convex 
function r/, G C^(M;M) by a sequence of functions of the form: r]l"\-) = YliPi"\' ~ ^1"^)^: 
with > 0. 

Entropy solution was first introduced by Kruikov [22] as the only physically admissible 
solution among all weak (distributional) solutions to scalar conservation laws. These weak 
solutions lack the fact of being unique for it is easy to construct multiple weak solutions to 
Cauchy problems (|2.2p . see |25j. 

Our next definition concerns classical sub-/super-solution to scalar conservation laws. This 
kind of solutions are shown to be entropy solutions, for the details see lemma [3?3l 

Definition 2.14 (Classical solution to scalar conservation laws) 

Let F{x,t,v) = g{x,t)f{v) with g, g^ G Lf^^{QT;R) and f G C^(R;M). A function v G 
W^''^{Qt) is said to be a classical sub-solution of (KM) "^^^^ v^{x) = t>(x,0) if it satisfies 

vt{x,t) + {F{x,t,v{x,t)))^ <0 a.e. in Qt- (2.19) 

Classical super-solutions are defined by replacing "< " with "> " in liS.19\) . and classical solutions 
are defined to be both classical sub- and super- solutions. 
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We move now to some results on entropy solutions depicted from [22]. 

Theorem 2.15 (Kruzkov's Existence Theorem) 

Let F, be given by Definition \2.1I\ and the following conditions hold: 

(EO) /GCi(M), 
(El) 5, G Cfe(QT), 
(E2) G C{Qt), 

then there exists an entropy solution v E L°°{Qt) of 112.^) . 

In fact, Kruikov's conditions for existence were given for a general flux function [22j Section 
4]. However, in Subsection 5.4 of the same paper, a weak version of these conditions, that 
can be easily checked in the case F{x,t,v) = g{x,t)f{v) and (E0)-(E1)-(E2), is presented. 
Furthermore, uniqueness follows from the following comparison principle. 

Theorem 2.16 (Comparison Principle) 

Let F be given by Definition \2. 12\ with f satisfying (EO), and g satisfies, 
(E3) g G W^'^iQr). 

Let u{x,t), v{x,t) € L°°{Qt) be two entropy sub- /super- solutions of Ii2.2\} with initial data 
G L°°(M). Suppose that, 

u°{x) < v°{x) a.e. m R, 

then 

u{x,t) < v{x,t) a.e. in Qt- 
Proof. See Section [U Appendix. □ 

It is worth noticing that in [22], the proof of the existence of entropy solutions of (|2.2p is 
made through a parabolic regularization of (|2.2p and passing to the limit, with respect to the 
L^ convergence on compacts, in a convenient space. 

At this stage, we are ready to present a relation that sometimes hold between scalar con- 
servation laws and Hamilton-Jacobi equations in one-dimensional space. 

2.3 Entropy- Viscosity relation 

Formally, by differentiating (12. ip with respect to x and defining v = Ux, we see that (|2.ip is 
equivalent to the scalar conservation law (|2.2I) with = and the same F. This equivalence 
of the two problems has been exploited in order to translate some numerical methods for 
hyperbolic conservation laws to methods for Hamilton-Jacobi equations. Moreover, several 
proofs were given in the one dimensional case. The usual proof of this relation depends strongly 
on the known results about existence and uniqueness of the solutions of the two problems 
together with the convergence of the viscosity method (see [8l[23l[2^). Another proof of this 
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relation could be found in [5] via the definion of viscosity/entropy inequalities, while a direct 
proof could also be found in ^2ij using the front tracking method. The case of a Hamiltonian 
of the form (12.31) is also treated even when g{x, t) is allowed to be discontinuous in the (x, t) 
plane along a finite number of (possibly intersected) curves, see |29| . 

In our work, the above stated relation will be successfully used to get some gradient es- 
timates of K. Although several approaches were given to establish this connection, we will 
present for the reader's convenience, a proof similar to that given in [8i Theorem 2.2]. For 
every Hamiltionian/flux function F = gf and every £ Lip{M.), let 



£V = {(VO), (VI), (V2), (EO), (El), (E2), (E3)}, 



in other words, 



The set of all conditions on / and g ensuring the 
existence and uniqueness of a Lipschitz continuous viscosity 
solution u S Lip{QT) of (|2.ip . and of an entropy 
solution V G L°°{Qt) of with = ul £ L°°(M). 



Theorem 2.17 (A link between viscosity and entropy solutions) 

Let F = gf with g e C^{Qt), G Lip{R) and £V satisfied. Then, 

V = Ux a.e. in Qt- 

Sketch of the proof. Let e > and 6 > 0. We start the prove by making a parabolic 
regularization of equation (12. Ih and a smooth regularization of uq and we solve the following 
parabolic equation: 

iuf + F{x,t,nf) = eu%i in Mx(0,r), 
\u^'\x,0) =u^^^ix) in M. 

For the sake of simplicity, we will denote u"^'^ by w and u^'^ by . Note that the first equation 
of (|2.20l) can be viewed as the heat equation with a source term F. Thus, we have: 

F\w](x,t) in Qt, 
^ ^ (2.21) 
in M, 

with F['w]{x,t) = F{x,t,Wx{x,t)). From the classical theory of heat equations, since F[w] £ 
^oc(Qt) and w° £ W^/oc (1^), there exists a unique solution w of (|2.21|) such that 

w e W^p'^(J^) CC Qt and l<p< oo. 

2 1 

Here the space Wp' (0), p > 1 is the Banach space consisting of all functions w € LP{Q,) 




having generalized derivatives of the form wt and Wxx in L^{^). For more details, see [24l 
Theorem 9.1]. We also notice that the space Wp'^{Q,) is continuously injected in the Holder 
space C"'"/^(r2) for a = 2 — I and p > |, see [25 . We use now a bootstrap argument to 
increase the regularity of w, taking in each stage, the new regularity of Flw] and the regularity 
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of w^. Finally, we get that w G C'^'^(M x [0, T)) (three times continuously differentiable in 
space and one time continuously differentiable in time). Prom the maximum principle and the 
L^'-estimates of the heat equation, see [MIE], it follows the uniform bound of u^''^ in W^^^{Qt), 
for p > 2. Therefore, we get as 5 ^ and e — > that: 

u^'^^u in C(Mx[0,r)), 

with u{x,0) = . We now make use of the stability theorem, [2l Theoreme 2.3], twice on the 
equation (|2.20l) to get that the limit u is the unique viscosity solution of (|2.ip . Hence, we have 
for any G C^{Qt) 

lim / / u'L'^chdxdt = — lim / / u'^'^Sx dx dt 

/ u(f)xdxdt= I I Ux4>dxdt. 
JR Jo Jr 

The appearance of Ux follows since u G Lip{QT). Moreover, as a regular solution, the function 
ye,s _ ^e,"? ggjygg ^]^g derived problem 



v'/ + (Fix, t, v''^))x = evli in R X (0, T), 
v^'^^x^Q) =u^/{x) in M, 



(2.22) 



and, according to [22l Theorem 4], the sequence v^'^ converge in L\^^{Qt), as e — > and 5 — > 0, 
to the entropy solution v of (12. 2p . Then, for any (/> G C^{Qt), 



lim / / v'^'^cp dxdt = I I vcf) dx dt. 



/ / if^^6dxd.t = I 

Consequently, 

/ Ux(j)dxdt= / / v(f)dxdt, 
Jr Jo Jr 

and Ux = V a.e. in Qt- □ 



Remark 2.18 The converse of the previous theorem holds under certain assumptions (see f21\ 

W)- 

Remark 2.19 In the multidimensional case this one-to-one correspondence no longer exists, 
instead the gradient v = \/u satisfies formally a non-strict hyperbolic system of conservation 
laws (see /ITjfl^). 



Throughout Sections 3 and 4, p will always be the solution of the heat equation (jl.Sp . The 
properties of the solution of the heat equation with such a regular initial data will be frequently 
used, we refer the reader to [3j [13] for details. 
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3 The approximate problem 



In this section, we approximate (11.41) and we pose a more restrictive condition (see condition 
(11.8(1 ) on the gradient of the initial data than of the physicaly relevent one (|1.6p . We prove 
a result of existence and uniqueness of this approximate problem, namely Theorem 11.21 and 
the reader will notice at the end of this section that this restrictive condition is satisfied for all 
time, and this what cancels the approximation in the structure of (11.4(1 and returns it to its 
original one. Finaly we present the proof of Theorem 11.31 

For every a > 0, we build up an approximation function fa S C;[(IR) of the function ^ 
defined by: 



it X > a. 



faix) 



2a — X 



or + ar{x — ay 



otherwise. 



(3.1) 



Proposition 3.1 For any a > 0, let fa be defined by 113. 1\) and H G C^(M) be a scalar-valued 
function. If 

Faix,t,u) = -H{px{x,t))pxx{x,t)faiu) (3.2) 

and Hp G Lip{U.), then the Hamilton- Jacobi equation 

j Kt + Fa{x,t,K,^) = in Qt, 
\ k{x,0) = K°(x) in R, 

has a unique viscosity solution n G Lip{QT). 

Proof. The proof is easily concluded from Theorems 12.61 12.71 and Proposition 12.101 after 
checking that the conditions (V0)-(V1)-(V2) are satisfied with 

g{x,t) = -H{px{x,t))pxx{x,t). (3.4) 

The condition (VO) is trivial, while for (VI) , we just use the fact that H is bounded on compacts 
and the fact that \px{x,t)\ < | l/o^l Il°°(r) in Qt- For the condition (V2), the regularity of p and 
H permits to compute the spatial derivative of g in Qt, thus we have: 

gx = -{H {px)plx +H{px)pxxx)- 

The uniform bound of the spatial derivatives, up to the third order, of the solution of the heat 
equation, and the boundedness of H on compacts gives immediately (V2). □ 



In the following proposition, we show a lower-bound estimate for the gradient of k obtained 
in Proposition 13.11 It is worth mentioning that a result of lower-bound gradient estimates 
for first-order Hamilton- Jacobi equations could be found in [Ml Theorem 4.2]. However, this 
result holds for Hamiltonians F{x, t, u) that are convex in the w- variable, using only the viscosity 
theory techniques. This is not the case here, and in order to obtain our lower-bound estimates, 
we need to use the viscosity /entropy theory techniques. In particular, we have the following: 

Proposition 3.2 Let G G C^(M;M) satisfying the following conditions: 
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(Gl) G{x) > G(0) > 0, 
(G2) G" > 0. 
Moreover, let 

H = GG' and < a < G(0). 

// satisfies: 

^^l{x) > G(/9°(x)), a.e. in M, 
then the solution k obtained from Proposition \3. 1\ satisfies: 

Kx{x,t) > G{p^{x,t)) a.e. in Qt- (3.5) 

In order to prove Proposition 13.21 we first show that G{px) is an entropy sub-solution of 

[ + (F(x,t,u;))^ = in Qt, , ^ 

< „ (3.6) 

\uj{x,{)) =uj^{x) in M, 

with = G{p^) and F is the same as in (13.2(1 . Before going further, we will pause to prove a 
lemma which makes it easier to reach our goal. 

Lemma 3.3 (Classical sub-solutions are entropy sub-solutions) 

Letv S W^'°°{Qt) he a classical sub-solution of (KM) ^^^^ v^{x) = v{x,0), then v is an entropy 
sub-solution. 

Proof. Let r^j, <I>, h and (p be given by Definition 12.121 Multiplying inequality (I2.19P by vlj^{v)4> 
does not change its sign. Hence, after developing, we have: 

v'i{v)vt(l) + r]'i{v)gxf{v)(l) + 7]'i{v)gf' {v)vx(l} <0, a.e. in Qt, (3.7) 

and since v is Lipschitz continuous, we use the chain-rule formula together with (|2.17p to rewrite 
(ITTll as: 

{7]i{v))t(l) + gJ{v)r][{v)(l) + gmv))x(l)<0, a.e. in Qt- (3.8) 

Upon integrating (13.8(1 over Qt and transferring derivatives with respect to t and x to the test 
function, we obtain: 



r]i{v{x,t))(j)t{x,t) -\- ^{v{x,t))g{x,t)(j)x{x,t)-\- 

h{v{x,t))gx{x,t)(l>{x,t)]dxdt+ / r]i{v°{x))(j){x,0)dx > 0, (3.9) 
which ends the proof. □ 

Following same arguments, classical super-solutions are shown to entropy super-solutions. We 
return now to the function G{px) and we are ready to show that it is indeed an entropy sub- 
solution of ((3.6p . In particular, we have the following: 



Lemma 3.4 The function G{px) defined on Qt is a classical sub-solution of h3. 0^) with initial 
data G{p^), hence an entropy sub-solution. 
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Proof of Lemma 13. 4L First, it is easily seen that G{pcc) £ W^'°°{Qt). Define the scalar 
valued quantity B on Qt by: 

B{x, t) = dt{G{p^{x, t))) + d^{F{x, t, Gipjx, t)))). 

Since < a < G(0), we use (Gl) to get fa{G{px)) = 1/G{px) and we observe that, 

B = G {px)pxt - d^. ' 

= G {px)Pxxx ~ 



G{px) 

G{px)[H' {px)plx + H{px)pxxx] - {G'{px)plxH{px)) 



\ G^ipx) J 

G{px)pxxx{G{px)G'{px) - H{px)) - pUH'{px)G{px) - H{px)G'{px)) 

GHPx) 

= -pIxG"{px)- 

The condition (G2) gives immediately that B <{). This proves that G{px) is a classical sub- 
solution of equation (13.6(1 and hence an entropy sub-solution. □ 



Proof of Proposition 13.21 From the definition of H and the properties of p, it is easy to 
check that g G C'^{Qt) and that £V is fully satisfied. Hence, we are in the framework of 
Theorem 12. 171 with = kP. This theorem gives that Kx is the unique entropy solution of (13. 6p 
with = K^. Moreover, by the previous lemma, G{px) is an entropy sub-solution of (13. 6p . 
Since 

i4 > G{pI), a.e. in M, 



we can apply the Comparison Theorem 12.161 to get the desired result. □ 

It is worth notable here that we do not know how to obtain the lower-bound on the spatial 
gradient Kx using the viscosity framework directly. However, for the case of the upper-bound, 
we can do so (see Remark l4.ip . At this stage, fix some e > 0, and let 



G%x) = Vx2 + e2 and a = G^(0) = e. 
It is clear that G'^{x) satisfies the conditions (G1)-(G2) with 

H\x) = X, 

and the Hamiltonian F from (13.21) takes now the following shape: 

F^{x,t,u) = -px{x,t)pxxix,t)fe{u). (3.10) 

Moreover, we have the following corollary which is is an immediate consequence of Propositions 
O and [321 

Corollary 3.5 There exists a unique viscosity solution n G Lip{QT) of 

( Kt + F^{x,t,Kx) = in Qt, 
\ k{x, 0) = k° G Lip{R) in R, 
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with satisfies: 



4 > ViPx? + a.e. in M. (3.12) 
Moreover, this solution k satisfies: 



Kx > V Px + C'-s- Qt- (3.13) 
The following lamma will be used in the proof of Theorem 11.21 

Lemma 3.6 Let c be an arbitrary real constant and take ip E Lip(M; M) satisfying: 

i^x^c a.e. in M. 

IfC^ C^(M;M) is such that ip — C has a local maximum or local minimum at some point xq S M, 
then 

Cx{xo) > C. 

Proof. Suppose that tp — ( has a local minimum at the point xq] this ensures the existence of 
a certain r > such that 

{ip - C){x) > - C){xo) \/x; \x - xo\ < r. 

We argue by contradiction. Assuming Cx{xo) < c leads, from the continuity of (x, to the 
existence of r G (0, r) such that 

Cx{x) < c Vx; \x — xqI < r . (3-14) 
Let yo be a point such that \yo — xo\ < r and yo < xq. Reexpressing (|3.14p . we get 

(C-cx)^(x)<0 Vx G (yo,a;o), 

and hence 

[{ip - cx)x{x) - (C - cx)xix)]dx > 0, 



'yo 

which implies that 

(^-C)(xo) > (V-C)(yo), 

and hence a contradiction. We remark that the case of a local maximum can be treated in a 
similar way. □ 

Now, we are ready to present the proofs of the first two theorems announced in Section [H 

Proof of Theorem [US Let K £ Lip(Qx) be the solution of (13. lip obtained in Corollary 
13.51 Let us show that it is the unique viscosity solution of (II. 4p among those verifying (|3.13p . 
To do this, we consider a test function (j) G C^{Qt) such that k — (j) has a local minimum at 
some point (xo,to) S Qt- Proposition 12.101 together with inequality (|3.13l) gives that 

K(.,to) G Lip(M) and i^x{-,to) >e a.e. in M. 

We make use of Lemma \3M with ■0(-) = ^(-i^o) and ("(•) = </'(•) ^o) to get 

4>xixo,to) > e. (3.15) 
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Since k is a viscosity super-solution of 

Kt - ft{K,x)PxPxx =0 in Qt, 

we have 

(l>t{Xf),tQ) - fei4'x{xo,to))px{xo,to)pxxixo,to) > 0. 

However, from (|3.15|1 . we get 

(t>t{xo,to)(l>x{xo,to) - Px{xo,to)pxx{xo,to) > 0, 

and hence k is a viscosity super-solution of 

KtKx = PxPxx in Qt- 

In the same way, we can show that k is a viscosity sub-solution of the above equation and hence 
a viscosity solution. The uniqueness of this solution comes from the uniqueness of the viscosity 
solution of (|3.11l) by reversing the above reasoning. □ 

Remark 3.7 Notice that the first equation of il-4\ ) can be viewed as a Hamilton- Jacobi equation 
of the type 

F{X,Vk)=0 in Qt, 
where F : Qt x i-^ M defined by: 

FiX,p)=pip2-pxiX)p,AX), 
with X = (x, t) and p = ipi,P2)- 

Proof of Theorem 11.31 Let 9 = Kx- By Theorem 12.171 is the unique entropy solution of 

Ot = {pxPxxfe{0))x in Qt, 
e{x,Q)=e^{x) in M, 

with 

e^{x) = kI{x) > V(pg)2 + e2, a.e. in R. 
Moreover, from Corollary 13.51 we have 

6 > A/pl + e^ a.e. in Qt, 

from which we deduce that fe{0) = ^ and hence our theorem holds. □ 

4 Proof of Theorem 11.61 

We turn our attention now to Theorem 11.61 Let < e < 1 be a fixed constant and take 

K^^'{X) = k\x) + €X. (4.1) 
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It is easy to check that the function k^'*^ belongs to Lip(M), and by condition p.6|) we get for 
a.e. X G R, 



Prom Theorem 11.21 there exists a family of viscosity solutions £ Lip^Qx) to the initial value 
problem (II. 4p that satisfy: 

— VPx + ^ a-e- in Qt- 

We will try to extract a subsequence of that converges, in a suitable space, to the desired 
solution 

4.1 Gradient estimates. 

Uniform bounds for the space-time gradients of will play an essential role in the determina- 
tion of our subsequence. 

I. e-uniform upper-bound for nf. 



Starting with the time gradient, we have for a.e. {x,t) G Qt- 

Kl{x,t)K%{x,t) = p^{x,t)p-^^{x,t), (4.2) 

and 

Kl{x,t) > ^/pl{x,t) + e^ > a.e. in Qt- (4.3) 

If pxix,t) = for some Lebesgue point {x,t) of and k^, it follows from (|4.2|) and (|4.3p that 
Kl{x,t) = 0. Otherwise, and since by (14. 3|) k| > \px\, we conclude that: 

\i^t\ < \\Pxx\\l°°{r) a.e. in Qt, (4.4) 

and hence we obtain an e-uniform bound of nf. 

For the space gradient, we argue in a slightly different way. The key point for obtaining 
the uniform bound of kI was the minoration of k% by \px\ so, roughly speaking, if we want to 
follow the same previous steps using the symmetry of (|4.2p in and one should also have 
an appropriate minoration of |«;^| by a well controlled function which no longer exists. 

II. Formal calculus and best candidate. 

We seek to find the best candidate to be an upper-bound of For this reason, we regard 
formally what is happening at the maximum of Dividing both sides of (|4.2|) by k| and 
differentiating with respect to the spatial variable, we get: 

e Pxx ~^ PxPxxx l^xxPxPxx , . rx 

i^x \'^x) 

Notice that k%x = at the maximum of k%. Multiplying equality (|4.5p by k% and integrating 
between and t, we obtain: 

d 



Jo It (^*^'^^''^) " / ^'^^"^ Pxpxxx)dT, 
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then ^ 

Jo 

and hence, 

\kI\ < \/2cit + C2, 

where 

ci = II(^'L)^IIl°°(m) + IIpSIIl°°(m)II/'LxIIl°°(r)) 

and 

C2 = (||k2||loo(r) + 1)2. 

The reason of taking C2 as above easily follows since kS'*^ = + e, by taking e small enough, 
namely less than 1. 

III. e-uniform upper-bound for 

Define the function S by: 



S{x,t) = \/2c\t + C2. 

Let us show that 5 is an entropy super-solution of (13. 6p with F given by (|3.10l) and uP{x) = 
S'(a;,0). Indeed, it remark that S E PF^'°°((3t), and we know that for every (x, t) G we 
have, 

5(x,i) > ^ = ||K°||ioc(R) + 1 > e, 

then 

/.(5(a;,t)) = ^^^ V(x,t)GQT. (4.6) 



The regularity of the function S permits to inject it directly into the first equation of (13.6 
Therefore, using (|4.6p . we have 



S, 



PxPxx \ Cl p1^ + PxPxxa 



S ) X \J2c\t + C2 V^Clt + C2 

^ Cl - {pix + faPxa'aO 
a/2C]T+~C2 

> 0, 

which proves, by Lemma [331 that S is an entropy super-solution of (13. 6p . Prom the discussion 
of the proof of Proposition 13.21 we know that k| is an entropy solution of (13.6(1 hence an entropy 
sub-solution. Since for e < 1 and a.e. x S M, we have. 



'^x (^) — 

< VS^ = 5(x,0), 



then we can use the Comparison Theorem 12.161 of scalar conservation laws to obtain: 



k\.{x, i) < Vcit + C2 < \/ c\T + C2 a.e. in Qt, (4.7) 
and hence we get an e-uniform bound for 



22 



Remark 4.1 We were able to obtain this e-uniform upper-bound of k% by using the viscosity 
theory techniques. In fact, we claim that C^''^{x,y,t) = K'^{x,t) — K^{y,t) and C'^{x,y,t) = 
{x — y)S{t) are two viscosity sub-/super-solutions of the following Hamilton- J acobi equation: 

dw 

= F{x,t,Wx) - F{y,t,-Wy) in V = {{x,y,t); x > y and t > 0} 

with initial data ('^''^{x, y, 0) = — and ("^{x, y, 0) = (x — y)S{0) respectively. Here 

F is given by 113.1 0\) . The claim is easy for C'^, and we refer to fSi when k'^ is a continuous 
viscosity solution of 113. We also notice that: C^''^{x,y,0) < C'^{x,y,0) y{x,y,0) £ V, and 
C^''^{x,y,t) = ('^{x,y,t) = for x = y, t > 0. Moreover, since C^'*^ and are continuous 
functions, we use the comparison principle of viscosity solutions (see J^) to obtain: 

k'{x, t) - K'{y, t)<{x- y)S{t) V(x, y, t) G P, 

hence, the estimate ( (^. ?| j holds. 

4.2 Local boundedness in VF^'°°. 

We now show that the family (K'^)o<e<i is locally bounded in W^'°"{Qt)- Let KCC Qt be a 
compactly contained subset of Qt, and {x,t) S K. Since k"^ is Lipschitz continuous, we can 
write, 

\K%X,t)-K^'^m<Chp\{x,t)\, 

where Cf-^ is the Lipschitz constant of which is independent of e from the previous estimates, 
namely (Ol) and (f477ll . Call this constant C. From the definition of k°'^(0) given by ([4T]l . it 
follows that, 

\K'{x,t)\ < C\ix,t)\ + \K\0)\, 

< C max |(2/,r)| + 1^0(0)1, 

{y,T)GK 

which is finite since K is bounded and hence, (K^)o<e<i is uniformly bounded in C{K). This, 
together with the uniform gradient estimates, gives the local boundedness of Kf^ in W^''^ {Qt) ■ 



4.3 Proof of theorem 11.61 

At this point, we have the necessary tools to give the proof of Theorem 11.61 We first recall 
that is a viscosity solution of an equation of the type (|4.2p . with a Hamiltonian independent 
of e (see Remark 13.71) and kP''' — > locally uniformly in M. By Ascoli's Theorem, there is 
a subsequence, called again Kf- , that converges to ac S Lip{QT) locally uniformly, and by the 
stability theorem (see [21, Theorem 2.3]), k is a viscosity solution of the initial value problem 

K-ti^x = PxPxx m. Qt, , , 

n (4.8) 
k(x,0) = K°(j;) in M. 

To end the proof, we still have to show the inequality 

>^x ^ \Px\ a-e. in Qt- 
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Again by Theorem 11.21 our Kf^ verifies for a.e. G Qt, 



> \Px(x,t)\, 

then for {y,t), {x,t) € Qt close enough, with a continuous function, the following inequality 
hold 

.%,t)-.-(x,t) ^ 
x-y 

Using the local uniform convergence of to k, we get a similar inequality with replaced 
with K and hence 

i^x ^ \Px\ a-c. in Qt- 

□ 



5 Problem with boundary conditions 

In this part of the paper, we deal with the same problem structure but with boundary conditions 
of the Dirichlet type. This sort of boundary conditions arises naturally in a special model of 
dislocation dynamics and will be explained in the following subsection. Our notations are kept 
untouched; the terms 6~^, 0~ , p and k still have the same physical meaning, while the domain 
is changed into the open and bounded interval 

/ = (o,i), 

of the real line. Although this problem seems to be an independent one, we will try to benefit 
the results of the previous sections by considering a trick of extension and restriction, in order 
to apply some of the previous results of the whole space problem. 



5.1 Brief physical motivation 

To illustrate some physical motivations of the boundary value problem, we consider a con- 
strained channel deforming in simple shear (see [E]). A channel of width 1 in the x-direction 
and infinite extension in the y-direction is bounded by walls that are impenetrable for disloca- 
tions (see Figure 1). The motion of the positive and negative dislocations corresponds to the 
x-direction. This is a simplified version of a system studied by Van der Giessen and coworkers 
[6], where the simplifications stem from the fact that: 

• only a single slip system is assumed to be active, such that reactions between dislocations 
of different type need not be considered; 

• the boundary conditions reduce to "no flux" conditions for the dislocation fluxes at the 
boundary walls. 

The mathematical formulation of this model, as expressed in [TS], is the system (|1.2p posed on 
Ix (0,r): 



a..^(M)-a.(nM)(' g;:;;j;^-:;:;;[ ))^o, 

Oti^,t)-Oxi^,t) 

9+{x,t) +9-{x,t) 



9,0-(x,t)+a.(e-(x,t)(^^^^l|— EM ) ) =0. 



(5.1) 
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Figure 1: Geometry of a constrained channel 

To formulate lieuristically the boundary conditions at the walls located at x = and a; = 1, we 
note that the dislocation fluxes at the walls must be zero, which requires that 



d^{e+-e-) = 0, at xG{0,l}. (5.2) 
Rewriting system (15. ip in a special integrated form in terms of p, k and we get 

1p. = *. 

Using (15. 2|) into the system (15. 3|) . we can formally deduce that p and k are constants along the 
boundary walls. Therefore, the remaining of this paper focuses attention on the study of the 
following coupled Dirichlet boundary problems: 

Pt = Pxx, in I X (0, oo), 

p{x,0)=p\x), in /, (5.4) 

p{0,t) = p{l,t) =0, VtE[0,cx)), 

and 

KtKx = PtPx, in I X (0,r), 

k{x,0) = K°{x), in /, (5.5) 

^ k(0, t) = k{0, 0) and k(1, t) = k(1, 0), Vt G [0, T]. 

Denote It by: 

/t = /x (o,r). 

There are two natural assumptions concerning p^ and k^, the first one is again the positivity of 
the dislocation densities 9^ and 6~ at the initial time, which yields to the following condition: 

> \P% (5.6) 



25 



and the second one has to do with the balance of the physical model that starts with the same 
number of positive and negative dislocations. In other words, if n"*" and n~ are the total number 
of positive and negative dislocations respectively at t = then: 

pO(l)-pO(0) = f pl{x)dx, 
Jo 

= [ {e+{x,o)-e'{x,o))dx, 

Jo 

= — n~ = 0, 

this shows that p^{l) = p^{0) and this is what appears in (15.41) . Up to now, formal relations 
between the initial conditions are only expressed. Whereas, required regularity, together with 
the announcement of the main results will be stated in the next subsection. 



5.2 Statement of the main results on a bounded interval 

Prom now on, the reader should not be confused with the term p that will always be the unique 
solution of the classical heat equation (|5.4p . The two main theorems that we are going to prove 
are: 

Theorem 5.1 (Existence and uniqueness of a viscosity solution) 

Let r > and e > be two constants. Take kP £ Lip{I) and p^ £ Cq°{I) satisfying: 

> G{pl) a.e. in I, 

where 

G{x) = + e2, 

then there exists a viscosity solution n G Lip{lT) of 15. 5\) . unique among those satisfying: 

i^x ^ G{px) a.e. in It- (5-7) 
Theorem 5.2 (Existence of a viscosity solution) 

Let T > and hP G Lip{I). Under the condition 115. 6)) satisfied a.e. in I, there exists a viscosity 
solution K G Lip{lT) of S5. 5\) satisfying: 

Kx ^ \Px\, fl-G- in It- 



5.3 Preliminary results 

Before proceeding with the proof of our theorems, we have to introduce some essential tools 
that are the core of the "extension and restriction" method that we are going to use. 



Extension of p over M x [0,T]. 

Consider the function p defined on [0, 2] x [0, T] by 



r'"-;' (-'O^'-. ,5.8) 

I — p{2 — X, t) otherwise. 
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this is just a antisymmetry of p with respect to the line x = 1. The continuation of p 
to ffi X [0, T] is made by spatial periodicity of period 2. A simple computation yields, for 
{x,t) G (1,2) X (0,r): 

pt{x,t) = -pt{2 - x,t) and pxx{x,t) = -pj;x{2 - x,t), 

and hence it is easy to verify that p |[i,2]x[o,T] solves (15.41) with / replaced with the interval (1, 2) 
and p^ replaced with its symmetry with respect to the point x = 1; the boundary conditions 
are unchanged and the regularity of the initial condition is conserved. To be more precise, we 
write down some useful properties of p. 

Regularity properties of p. 

Let r and s are two positive integers such that s < 2. From the construction of p and the above 
discussion, we get the following: 

i) pt and px are in C(IR x [0,T]), 

ii) rt = Oon Zx [0,ri, 

' ^ ^ (5.9) 

iii) /5t = /5a.x on (M\Z) X (0,T), 

iv) ||9r5^p(-,0llL-{M)<C, VtG[0,r], 

Where C is a certain constant and the limitation s < 2 comes from the spatial antisymmetry. 
These conditions are valid thanks to the way of construction of the function p and to the max- 
imum principle of the solution of the heat equation on bounded domains (see [3l [13]). 

Let 

g{x,t) = -pt{x,t)pxix,t). (5.10) 

Prom the above discussion, it is worth noticing that this function is a Lipschitz continuous 
function in the x-variable. 

The following three lemmas will be used in the proof of Theorem 15.11 

Lemma 5.3 (Entropy sub-solution) 

The function G{px) is an entropy sub-solution of 



Wt + {gfe{w))x = 0, in Qt, 
w{x,Q) = ViP {x) in M, 



(5.11) 



where fe is given by hS. and w^{x) = G{px{x,0)). 

Proof. Similar to Lemma \3M □ 
Lemma 5.4 (Differentiability property) 

Let u{x,t) be a differentiable function with respect to {x,t) a.e. in Qt- Define the set M by: 

M = {x G M; u is differentiable a.e. in {x} x (0,T)} , 
then M is dense in M. 
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Proof. Define n G N to be the Lebesgue n- dimensional measure. Let N C Qt be the set 
defined by: 

N = {{x,t) G Qt] u is not differentiable on {x,t)} , 
and let In be the characteristic function of the set N. Since C'^{N) = 0, we can write, 

1n{x-, t)dxdt = 0. 

JQt 

Using Pubini's theorem we get 

cT 



I g{x)dx = 0, with g{x) = I / I 



N{x,t)dt > 



then 

g = a.e. in M 

and consequently 

J = {x; g{x) 0} verifies £^(J) = 0. 

In other words, 

Vx G M \ J, u{x, •) is differentiable with respect to a.e. in (0,T), 
hence M \ J C M which implies our lemma. □ 

In the next lemma, we show a lower-bound estimate for the gradient of k analogue to 
(|5.7|) . This was previously done for Kx in the case where 5 is a twice continuously differentiable 
function using mainly Theorems 12.171 and 12.161 Here, the way of extending the function p over 
Qt makes g loose some of the regularity stated in Theorem 12. 171 However, the following lemma 
shows that a similar result holds in the case g G W^'°°{Qt)- 

Lemma 5.5 The function G L°°{Qt) is an entropy solution of 115.11]) with initial data 
w° = kl£ L°°(M). 

Proof of Lemma 15.51 Let g be an extension of the function g on M? defined by: 

' g(x,t) if (a;,t)GQr, 
g{x,t) = I g{x,T) if t > T, (5.12) 
J{x,0) if t<0. 

Consider a sequence of mollifiers ^" in and let ^" = g*£,"- Remark that, from the standard 
properties of the mollifier sequence, we have g^ G C^{M?) and: 

(p ^ g uniformly on compacts in Qt, (5.13) 

and 

-g'i ^ g, in Lf^ ,(Qt), 1 < P < 00, (5.14) 
together with the following estimates: 

Wtdir\\L^iQr)<Wtdt,g\\L^^Q^^ for r,sGN,r + s<l. (5.15) 
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Now, take again the Hamilton-Jacobi equation (I5.27P with g replaced with g"^' 



ut + rfeiu^) = in Mx(0,r), 
n(x,0)=K°(x) in M, 



(5.16) 



and notice that the above properties of the function g" enters us into the framework of Theorem 
I2.17[ Thus, we have a unique viscosity solution k" € Lip(QT) of (|5.16ll with initial condition 
kP whose spatial derivative € L°^{Qt) is an entropy solution of the corresponding derived 
equation with initial data k^. Prom Remark 12.111 and (|5.15p . we deduce that the sequence 
(/*")n>i is locally uniformly bounded in W^'°^{Qt) and that: 



\'^x\\l°°(Qt) — \\'^x\\l°°{M.) +T\\gx\\L^(Qrp)\\fe\\L° 



(5.17) 



Moreover, from (I5.13p . we use again the Stability Theorem of viscosity solutions [21 Theorem 
2.3], and we obtain: 

— > K locally uniformly in Qt- (5.18) 
Back to the entropy solution, we write down the entropy inequality (see Definition l2.12p satisfied 



Qt 



r]{k'^)(Pt + HK'^W<Px + Kk]i)g^cl))dxdt+ / r]{kl)cj){x,0)dx>0, 



(5.19) 



where rj, h and cp are given by Definition 12.121 Taking (I5.17P into consideration, we use 
a property of bounded sequences in L^{Qt) (see [HI Proposition 3]) that guarantees the 
existence of a subsequence (call it again k") so that, for any function ip £ 



V'(k^) ^ C/^ weak-* in L°°{Qt). 
Purthermore, there exists /i G L°°{Qt x (0, 1)) such that: 

^{fi{x,t, a))da = U^{x,t), for a.e. {x,t) £ Qt- 



(5.20) 



(5.21) 



Applying (I5.20p with ^p replaced with rj, $ and h respectively, and using (I5.2ip . we get: 

viKi-)) ^ / il{fi{.,a))da weak-* in L°°{Qt), 
Jo 

Hf^xi-))^ I H^^{■,a))da weak-* in L°°{Qt), (5-22) 

JO 

KKi-))^ I h{ii{.,a))da weak-* in L°°{Qt). 
Jo 



This, together with (|5.13p . (|5.14l) permits to pass to the limit in (|5.19ll in the distributional 

sense, hence we get: 



It jo 



r]ip,i., a))(j)t + ^{K-^ ct))g(px + Kii{., a))gx(t)] dxdtda+ 



/ r]{kl)(pix,0)dx > 0. 
Jm. 



(5.23) 
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In [m Theorem 3], the function // satisfying (|5.23|) is called an entropy process solution. It 
has been proved to be unique and independent of a. Although this result in pj^j was for a 
divergence-free function g G C^(Qt), we remark that it can be adapted to the case of any 
function g G VF^'°°(Qt) (see for instance Remark 16.21 and the proof of [lH Theorem 3]). Using 
this, we infer the existence of a function z G L°°{Qt) such that: 

z{x,t) = fi{x,t,a), for a.e. {x,t,a) e Qt x (OA), (5-24) 

hence, z is an entropy solution of (IS.lip . We now make use of (I5.24p and we apply equality 
(|5.2ip for ^p{x) = x to obtain, 

z = weak-* lim in L°°{Qt)- (5.25) 

n— >oo 

Prom (|5.25l) and (I5.18P we deduce that, 

z{x,t) = kx{x,t) a.e. in Qt, 
which completes the proof of Lemma 15.51 □ 

5.4 Proofs of Theorems [SUl [O 

Proof of Theorem 15.11 We extend the function to G Lip(R) in the following way: 



r ,.0 



k^{x) 



K^ix) if XG[0, 1], 

/9°,||ioc(^) + e)(x- 1) + K°(1) if x>l, (5.26) 
A\\pI\\l--{i) + ^)x + k\0) if x<0. 



Consider the initial value problem defined by: 

(ut + gUux) = in Mx(0,r), 
\ u{x,0) = K°(x) in M. 

This is a Hamilton- Jacobi equation with a Hamiltonian F G C{Qt x M) defined by: 

F{x,t,u) = g{x,t)fe{u). 

Prom the regularity properties of p, we can directly see that (V0)-(V1)-(V2) are satisfied; 
this is quite similar to what was done in Proposition 13.11 Since is a Lipschitz continuous 
function, we deduce from Theorems 12.61 \2J\ and Proposition 12.101 the existence and uniqueness 
of a viscosity solution k G Lip{QT) of (I5.27p . Moreover, in order to recover the boundary 
conditions given by (15. 5p on dl x [0, T], we proceed as follows. Let M be the set defined by 
Lemma [EH and let x G M. For every t G [0, T], we write: 

\k{x,t)—k{x,0)\ < / \ks{x, s)\ds < / \F{x, s, kx{x, s))\ds < / {\F{0, s, kx{x, s))\ + C\x\) ds. 



In these inequalities we have used the fact that k is a Lipschitz continuous viscosity solution 
of (|5.27ll and hence it verifies the equation in Qt at the points where it is differentiable (see 
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for instance [2]). Also, we have used the condition (Fl) with p = q and Cr = C, a constant 
independant of R. Now from (|5.9l) -(ii). we deduce that: 

\F{0,s,kx{x,s))\ = \px{0,s)pt{0,s)fe{krc{x,s))\ = 0, for a.e. s G (0,t), 

and hence we get 

\k{x,t) - k(x,0)| < C\x\t. (5.28) 
Since M is a dense subset of M, we pass to the limit in (I5.28P as x — > and the equality 

k{0, t) = k{0, 0) = K°(0) Vt G [0, T] 

holds. Similarly, we can verify that k{l,t) = k{l,0) = K^{1) for all t G [0, T]. 

Existence. The extension k^ of outside the interval I is a linear extension of slope | l-L°°(/)+ 
therefore we have, 

f^'{-)>VMW+^=G{pl{-)), a.e. in M. (5.29) 

Prom Lemma 15.51 we know that kx is an entropy solution of equation (15.11(1 and from Lemma 
15.31 we know that G{px) is an entropy sub-solution of (I5.1ip . Since (I5.29P holds, we use the 
Comparison Theorem 12.161 to get, 

kx{x,t) > ^/pl{x,t) + e2 > e > 0, for a.e. {x,t) E Qt- (5.30) 

Take k to be the restriction of k on It where k^ and p have their automatic replacements 
and p respectively on this subdomain. It is clear that k £ Lip{lT) is a viscosity solution of: 

' K-t + giMx) =0 in It, 

< k(x,0) = K^{x) in /, (5.31) 

^ K(0,t) = K°(0) and K(l,t) = K°(l) V < t < T, 

where g{x,t) = —pt{x,t)px{x,t) and Kx{x,t) > G{px{x,t)) for a.e. {x,t) £ It- We also notice 
that K is a viscosity solution of (|5.5p . for it suffices to follow the same steps of the passage from 
the viscosity solution of (|3.1ip to the viscosity solution of (|1.4|) (see the proof of Theorem 11.21 
for details). 

Uniqueness. Since the function 

H{x,t,u) = gix,t)Mu) G CilT xR) 

satisfies for a fixed t: 

\Hix,t,u) - H{y,t,u)\ < C{\x - y\{l + |n|)), 

for every x, y £ (0, 1) and u E M, we use [21 Theorem 2.8] to show that k is the unique viscosity 
solution of (I5.3ip . We claim that k is the unique viscosity solution of (|5.5p . Indeed, we can 
also follow the same mechanism as in the proof of Theorem 11.21 □ 

We now move towards the proof of Theorem 15.21 that has the same flavor of what was done 
in Section [H We just need to care about the change in the structure of our problem and the 
boundary conditions. Our first step will be the following lemma. 
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Lemma 5.6 Let c\ and C2 be two positive constants defined respectively by: 
and 

Then the function S defined on Qt by: 

S{x,t) = \j2c\t + C2 

is an entropy super-solution of i5.ll]} with 

w^{x) = S{x,0) = \\k^\\l°°(i) + 1- 
Proof. See Subsection 5.1-III. □ 

Proof of Theorem 15. 2L Let e > be a fixed constant. Define kP''^ £ Lip{M.) by: 

^ K°{x)+ex if X G [0,1], 

(ll'^SllL-(/)+e)(^-l) + ('^°(l)+e) if x>l, (5.32) 
[{\\K-l\\Lo-{i)+e)x + k\0) if x<0. 



Since > \p^\ a.e. in /, it is clear that for a.e. x S M we have 



k'/ > Gipl), 

and hence, from the discussion of the proof of Theorem 15.11 there exists a unique viscosity 
solution S Lip^Qt) of 



^tf^l = PtPx in Qt 

k'{x,0) = G Lip{R) in M, 



oe, . - . ^ (5-33) 



unique among those satisfying: 

kl > G{px) a.e. in Qt- (5.34) 

Assume without loss of generality that e < 1. The e-uniform bound for kl is trivial, it suffices 
to use directly the equation satisfied by k'' together with (I5.34p . And the e-uniform bound for 



k% follows from Lemma 15. 61 and Theorem 12.161 since 

k|(x,0) < + e < + 1 = y/ci = S{x,Q). 

Following exactly the same technic of Section lU namely the proof of Theorem 11.61 we get that 
the sequence k"^ converges locally uniformly to k in Qt with k G Lip{QT) satisfies, 

kx ^ \Px\ Si.e. in Qt (5.35) 

and 

k{x,0) = ko{x) in M, (5.36) 
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where kq is the uniform limit of the sequence kP''^ in M. Theorem 15.11 guarantees that 

/i^(0,t) = K°'^(0) = kO(0), (5.37) 

and 

= = + e, (5.38) 

for all t E [0, T]. From (|5.37ll . (I5.38|) and the pointwise convergence, up to a subsequence, of 
to k, we deduce that 

k{0, t) = lim k'{0, t) = k°(0), Vt G [0, T], (5.39) 

and 

K(l,t) = limK"(l,t) = lim(K°(l) + e) = K°(l) VtG[0,T]. (5.40) 

Take k to be the restriction of k over It', p and have their automatic replacements p and k'' 
respectively on this restricted domain. Prom (|5.35ll . (I5.36p . (|5.39|) and (I5.40p . we deduce that 
K is the required solution. □ 

6 Appendix: Proof of Theorem 12.161 

We will work on the entropy inequality (12.18(1 satisfied by u and its analogue satisfied by v, 
using the dedoubling variable technique of Kruzhkov (see p2|) and following the same steps 
of [m Theorem 3], taking into consideration the new modifications arising from the fact that 
we are dealing with sub-/super-entropy solutions and the fact that g € W^'°^{Qt) is not a 
gradient-free function. 

The proof can be divided into three steps. Denote by Br = {x £ M; |x| < r} for any r > 0, 
F^{u, v) = sgn^{u - v){f{u) - f{v)), 

y~ = l|y||L-(QT) for every y € L°°(gT) (6.1) 

and 

Mf = max \f'{x)\. (6.2) 

x I <max(M°° ,1)°° ) 

In step 1, we prove that the initial conditions u^, satisfy for any o > 0: 



lim ^ /" [ {u{x,t) - vP{x))"^dxdt = 0, (6.3) 

'o JBa 



lim-/" I {v{x,t) -v'^{x))'dxdt = 0, (6.4) 
■^^OrJo JBa 

respectively. 

In step 2, The following relation between u and v is shown: 

[{u{x, t) - v{x, t))^iJt + F~^{u{x, t),v{x, t))g{x, t)^'^] dxdt > 0, (6.5) 
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for every G Cq^M x (0,r);M+). 

After that, we define A{t) for < t < min(r, ^) and lo = g°^Mf, by: 



A{t)= {u{x,t) -vix,t))^ dx. (6.6) 



Ba 



In step 3, we show that A is non-increasing a.e. in (0, min(T, ^)) and we deduce that 

u{x,t) < v{x,t) a.e. in Qt- 



Step 1: Proof of dnSD, (iGTil) . 

Let be a sequence of molhfiers in M with = ^. Recall that the function ^ G Co°(]R) 
satisfies the following properties: 

supp(e) = {x G M, e(x)/0}cSi; 

e>0, e(-x)=C(x); 

(6.7) 



^{x)dx = 1; 

^"(x) = n^{nx). 
Let r G M such that < r < T and define the function 7 by: 

{T — t 
if < t < T, , , 
- - (6.8) 
if t > r. 

Take a > and a test function iIj £ C^(M;M+) such that, 

tp{x) = 1 for X G Ba- 

Let y G M be a Lebesgue point of and we make use of inequality (I2.18P with k = u^{y) 
and the test function (p{x,t) = ip{x)j{t)^'^{x — y) (this is possible since is a permissible test 
function). Integrating the resulting inequality with respect to y over M yields: 



with 



and 



Ti(n, r) + r2(n, r) + r3(n, r) + r4(n) > 0, (6.9) 

Ti{n, t) = -- [ I {u{x, t) - nO(y))+V(x)r (2; - y) dxdydt, (6.10) 
'T Jo Jr2 

T2{n,T)= [ I F+{u{x,t),u\y))g{x,t)-i{t){ij{x)C{x-y)).dxdydt, (6.11) 
Jo Jm? 

%{n,T) = - r [ sgn+{u{x,t) - u\y)) f {u^ {y)) 
Jo JR2 

5^(x, t)-f{t)^{x)e{x - y)dxdydt (6.12) 

r4(n) = / (u°(x) -u°(2/))+^(x)C(x-2/)dxdy. (6.13) 
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Using the change of variables: x = x ^ y = x — ^in (I6.10p . and denoting again by the 
new variables (x',y') yields: 

Ti{n,T) = -- [ [ ( (u{x,t)-v^(x-^y ^{x)i{y) dxdydt, (6.14) 
JO JBi JR ^ \ nJ J 

Using that, 

[u — v)'^ — [u — w)^ <{w — u)^ V«, w G M, (6.15) 



we infer that: 

r*(r) 



ri(n,r) + i / / {u{x,t)-u^{x))^i)(x)dxdt< 



T 



K,u JB 



n° - ^) - S.iy)dydx, (6.16) 



where is the support of ip. Same upper-bound, independent of r, could be obtained for 
T^{n). Furthermore, since G L°°(R), thus integrable over K^, we use the Lebesgue differen- 
tiation Theorem to show that the right side of (|6.16|) tends to when n becomes large. Now, 
let e > 0, 3no such that 

Ti(no,r) +r*(T) < I and nino) < |, Vr > 0. (6.17) 

We also remark that the integrands of the right hand sides of (16. lip and (I6.12P are bounded 
and hence, for this particular no we can choose some tq such that V < r < tq, we have: 

T2{no,T)<^ and T3(no,T)<|. (6.18) 

Prom (fHTTD . (I6T8II and ([HJl), we infer that, 

0<T*(r)<e, VO<r<ro. 

Since ipi^x) = 1 over Ba, (|6.3p is proven. Arguing in the same way, we can prove (16. 4p . The 
slight difference is using a similar inequality of (I6.15P with (•)"'" replaced with (•)". 

Step 2: Proof of i^^. 

It suffices to prove (|6.5p for any function ip £ C^{Qt',^+)- We may also assume, without 
loss of generality, that there is some c > such that tpixjt) = for t S (0, c) U (T — c, T). For 
n > i, let be the usual mollifier sequence in M and consider the function (f>{x, t, y, s) defined 
for {x, t) e Qt and {y, s) G Qt by, 

^{x, t, y,s) = i; ^-4A C{x - y)Cit - s). 



The function (p hence satisfies 

<Pi.,.,y,s)eC^{QT;^+) and <p{x,t, ., .) e C^iQT;^+). 
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Fix some (y, s) G Qt for which the function v is well defined (this is valid almost everywhere). 
Since u is an entropy sub-solution of (12. 2p . we consider the relation (I2.18|) satisfied by u with 
k = v{y,s) and the test function (p{., .,y, s). Upon integrating this inequality with respect to 
{y, s) over Qt, we get: 



{{u{x,t) - v{y,s))^(/)t{x,t,y,s) + {u{x,t),v(y, s))g{x,t)(l)x{x,t,y, s) 
-sgn^{u{x,t) - v{y,s))f{v{y,s))gxix,t)(j){x,t,y,s)} dxdtdyds > 0. (6.19) 



Similar inequality could be obtained since v is an entropy super-solution of (12. 2p . We just swap 
+ , u and {x,t) with — , v and {y,s) respectively, hence: 



{{v{y,s) -u{x,t)) (ps{x,t,y,s) + F {v{y,s),u{x,t))g(y,s)(j)y{x,t,y,s) 



-sgn {v{y,s) - u{x,t))f{u{x,t))gxiy,s)(j){x,t,y,s)] dxdtdyds > 0. (6.20) 
Summing (I6.19P and (|6.2QI) and using the elementary identities: 

x~ = (—2;)^ and sgn~{x) = —sgn~^{—x), Vx G M, 

we get, for u = u{x,t) and v = v{y,s), 

Zi + Z2 + Z-i> 0, (6.21) 



with: 



Zi = {u — v)'^{4)t + 4's)ix,y, t,s) dxdtdyds, (6.22) 
F^{u, v) [g{x, t)(j)x {x, y, t, s) + g{y, s)(py{x, y, t, s)]dxdtdyds, (6.23) 



23= sgn+{u- v)[f{u)g^{y,s) - f{v)g^{x,t)]<j){x,y,t,s)dxdtdyds. (6.24) 



We now compute the first partial derivatives of the function (/>. For {x,t,y,s) S Qt x Qt, we 
have: 

Mx, t, y, s) =Cix- y) Qv't ^) e{t - s) 

+A^'-^)e\t-s)), (6.25) 



t, y, s) =e{x-y){ \i^t ^-^] eit - s) 



-i^[^^,'-^)cit-s)], (6.26) 



Mx, t, y, s) =e{t-s){ f ^ ) e{x - y) 
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(6-27) 



<t>y{x, t, y, .) = e{t - s) Qv'. ^ ) r (x - y) 

(6.28) 

Using these relations in (I6.2ip and performing the following change of variables, 

x =(x + y)/2, y = n{x -y), t = {t + s)/2, s = n{t - s); 

denote the new variables x' , t' , y' , s' by x, t, y, s and Q4 = Qt x Bf. Also, for the simplicity 
of expressions, denote 

This altogether yields: 



X = X -\ , = t -\ , X =x , t =t . 

2n 2n 2n 2n 



with: 



Xi + X2 + X3 + X4> 0, (6.29) 



Mi= {u{x^,t^) - v{x' ,t')yiJt{x,t)^{y)(,{s)dxdtdyds, (6.30) 

X2 = \l F+{u{x+,t+),v{x~,t-))ig{x^,t+)+g{x',t-))x 

'Q4 (6.31) 



2 



t/jxix, t)i{y)i{s)dxdtdyds, 

X3= [ F+{u{x+,t+),v{x-,t-)){g{x+,t+)-g{x~,t-))x 

Jq4 (6.32) 

ip{x,t)n^ {y)(,{s)dxdtdyds, 

Xi = f sgn+{u{x+,t+)-v{x-,t-))[f{u{x+,t+))gx{x-,t-)- 

Jq4 (6.33) 

f{v{x~,t~))gx{x+,t^)] ip{x, t)C{y)C{s)dxdtdyds. 
At this point, it is worth mentioning that we will frequently use the following Lemma from |23| . 

Lemma 6.1 IfTG Lip{M.) satisfies \T{u) — T{v)\ < Co\u — v\, then the function 

H{u, v) = sgn^{u — v){T{u) — T{v)) 

satisfies \H{u, v) — H{u ^v)\ < Co{\u — u'\ + \v — v'\ (see fSR Lemma 3]). 

Consider now (|6.30l) . Since {u — f )+ = sgn^{u — v){u — v), we make use of Lemma 16.11 to 
obtain: 

Xi — {u{x,t) — v{x,t))~^iljt{x,t)dxdt < 

Jqt 

|u(x+,t+) - u{x,t)\{ipt)°°^{y)i{s)dxdtdyds 

\v{x- ,t-) - v{x,t)\{iJt)^^{y)C{s)dxdtdyds \ , 
Ik^ JbI I 
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where, by the Lebesgue Differentiation/Dominated Theorems, the right hand side of this in- 
equality tends to as n — > oo, and hence: 

Xi ^ {u{x,t) — v{x,t))~^iljt{x,t)dxdt as n — > oo. (6.34) 

Jqt 

Let us now turn to (|6.31|) : using the fact that g G W^'°°{Qt) and hence Lipschitz continuous 
over the compact K^, and the fact that F~^{u, v) is Lipschitz continuous in u and v (see Lemma 
I6.ip . we get: 



X2 — {u{x,t),v{x,t))g{x,t)il}x{x,t)dxdt 



< 



g'^Mfilj^ { / \u{x+,t+) -u{x,t)\C{y)^{s)dxdtdyds 



(6.35) 



+ / / \v{x ,t ) — v{x,t)\S,{y)S,{s) > dxdtdyds 
Jk^ JbI J 

+ ^C7((5.)°°, {gtT, {i^ccr,Mf, u'^,v^,T), 

and also, by the Lebesgue Differentiation/Dominated Theorems, the left hand side of this 
inequality tends to as n ^ 00, hence: 

X2 ^ I F^{u{x,t),v{x,t))g{x,t)il)j.{x,t)dxdt as re — > 00. (6.36) 
Jqt 

We now study the two terms and X^. From the fact that g £ W^'^^iQr), we remark that 
for a.e. {x,t,y,s) G Qt x Qt, we have: 

g{x~,t') -5(2;+, t+) = gx{x~,t~){-y/n) + gt{x~ ,t~){-s/n) +o(- 

\n 

We also remark that the term gx{x~^,t^) in could be replaced with gx{x~,t~), since this 
adds a term that approaches as re becomes large. This term will be omitted throughout what 
follows and we denote the new by X^. From these two remarks, we rewrite X^ and X^ to 
get: 

X^= [ sgn+{u{x+,t+)-v{x-,t-)){f{u{x+,t+))-f{v{x-,t-))) 

jQi (6.37) 

{ygx{x',r) + sgt{x',t'))ip{x, t)^' {y)^{s)dx dt dy ds + £(n), 
where C{n) ^ as re ^ 00, and 



X^= sgn+iu{x+,t+)-v{x-,t-))ifiu{x^,t+))-fivix-,t-))) 

Iq^ (6.38) 

gx{x' ,t~)7p{x, t)^{y)S^{s)dx dt dy ds. 
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The term C{n) will also be omitted for simplification and we denote the new by X^. Let 



X^^ = X^ + X^, hence: 
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Q4 



F+(u(x+,t+),w(x ,t ))gx{x ,t )ip{x,t){y^{y)C{s))ydxdtdyds 



v2n 



(6.39) 



+ / F+(u(x+,t+),w(x ,t-)))gtix-',t~)i;{x,t){s^{y)^{s))ydxdtdyds. 

'Q4 

In A'g^l^ and X^^, the term ip{x,t) could be replaced with for this also adds a term 

getting small when n — > oo. We keep the same notations for X^^- Since y^{y)S,is) is 
a compactly supported smooth function in Q4, we have: 



S4 



F+(n(x ,t ,t ))gx{x ,t )'\l){x ,t ){yS,{y)^{s))ydx dt dy ds = 0. 



(6.40) 



Moreover, since F^{u,v) is Lipschitz continuous, we obtain: 



In 
34 



F+(u(a; ,t ),v{x ,t ))gx{x ,t )^{x ,t ){y^{y)i{s))ydx dt dy ds 



(6.41) 



'S4 

<Mf{gx)^^l^°° I / |'a(x+,t+) -n(x-,t-)|(ix(itd2/ds, 
where is the support of ip. Therefore, by the Lebesgue Differentiation/Dominated Theo- 



rems, we deduce that the right hand side of (I6.4ip tends to as n — > 00, hence we have: 



In 
34 



In a similar way we can show that 



X. 



2n 



as n — > cxD. 



as n — > 00. 



(6.42) 



(6.43) 

Prom (lOill . (lOHIl . (lOSjl and (fOHIl . passing to the limit in (lOnil yields (l63|), which con- 
cludes the proof of step 2. 



Step 3: u{x,t) < v{x,t) a.e. in Qt- 

Let us first show that the function A{t) defined in (|6.6p is non-increasing a.e. in (0, min(T, ^)). 
Take a > and recall that oj = g°°Mf; let < ti < t2 < min(T, ^), < e < 
min(ti, min(T, ^ — t2), and 6 > 0. Consider the function (/> G Co(M+,[0, 1]) such that 
(j){x) = 1 Vx e [0, a], 4>{x) = Vx G [a + 5, 00), and (f) < 0. Define by: 



r 



r,(t) 



i - (ti - e) 



1 

(^2 + e)-t 



if < t < ti - e 

if ti-e<t<ti 

if ti<t<t2 

if ^2 < t < ^2 + e 

if t2 + e < i < 00. 



(6.44) 
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One can take in (16. 5|) the permissible test function 

V'(x,i) = 0(|x| +iot)r^{t). 



This yields: 



Ei{5,e) 



{u{x, t) — v{x, t))~^(p{\x\ + ujt)dxdt ■ 

E2{5,e) 



I ft2+t 

{u{x, t) - v{x, t))^<t){\x\ + ijjt)dxdt > E{5, e), 



(6.45) 



with ^ 

E{6,e) =-/ [uj{u{x,t) -v{x,t))^ + sgn'^{{u{x,t) -v{x,t)))x 

Jo Jr (6.46) 

{f{u{x,t))-f{v{x,t)))^g{x,t)](t>'{\x\+uJt)r,{t)dxdt. 

\x\ 

We claim that E{d, e) > 0. Indeed, since < and > 0, it suffices to show that 

uj{u{x, t) — v{x, t))"*" + sgn'^{{u{x, t) — v{x, t))) x 

{fiuix,t))-f{v{x,t)))^gix,t)>0 a.e. in Qt- ^^'^^^ 

\x\ 

Two cases can be considered, either u{x,t) < v{x,t); in this case it is easy to verify (|6.47p . or 
u{x,t) > v{x,t); in this case we use, from the definition of a;, the fact that 

X 

{f{u{x, t)) - f{v{x, t)))-—rg{x, t) > -uj{u{x, t) - v{x, t)), 

\x\ 

hence our claim holds. Relation (|6.45l) now holds with E{5, e) replaced with 0. We regard the 
integrand term of Ei{6, e) in (I6.45P and we notice that for ti — e < t < ti, we have: 

{u{x, t) — v{x, t))~^(p{\x\ + LUt) = {u{x, t) — v{x, t))~^(p{\x\ + LOt)lAs, 

where I^/ is the characteristic function of the set Ag defined by: 

A's = {{x,t); ti - e < t < ti, <\x\+Lot < a + 5}. 
Remark that the set A'^ shrinks, as 6 becomes small, to 

a' = {{x, t); ti - e < t < ti, < |x| + cot < a} 
with (j){\x\ + Lot) = 1 over A. It is easy now to see that as 5 ^ 

{u{x,t) — v{x,t))^ (p(\x\ + Lut)!./ {u(x,t) — v{x,t))~^lA a.e. in Qt- 

s 

However, since {u{x,t) — v{x,t))^ £ L°°{Qt), we use the Lebesgue Dominated Theorem to get: 
Ei{5,€)^-[ I {u{x,t) -v{x,t))+dxdt as 6 ^0, (6.48) 
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in other words, 

rti 



1 

Ei{6,e)^- / A{t)dt as 6^0, (6.49) 



with A(t) given by (|6.6p . Similar arguments shows that: 

E2{S,e)^- A{t)dt as 5^0. (6.50) 

f Jt2 

Note that ^ G L^{0,T); let ti and t2 be Lebesgue points of the function A such that < ti < 
t2 < min(T, ^), one can easily deduce from (I6.49p . (I6.49P and (|6.45ll letting e tends to that 

A{ti) > A{t2), 

hence A is a.e. non-increasing. We use this property enjoyed by A to get the comparison 
principle. In fact, using the elementary identities: 

(m — w)^ < {u — w)'^ + {v — w)~ 
{u — v)^ < {u — w)~ + {v — w)'^ 

V n, f , ti; S M, we calculate for a.e. (x, t) S Qt '■ 

(n(x, t) - v{x, t))+ < {u{x, t) - u°(x))+ + {v{x, t) - v^{x))- + (n°(x) - t;°(x))+. 

Since vP{x) < v^{x) a.e. in R, we get for a.e. {x,t) G Qt'- 

{u{x,t) -v{x,t))^ < {u{x,t) - u°{x))+ + {v{x,t) - v°{x)y . (6.51) 

Using (I6.5ip . for r G (0, T), we calculate: 



- r A{t)dt <- r I 

T Jo T Jo JSa 



{u{x,t) — v{x,t))^ dxdt < 



-[ [ {u{x,t) - u^{x))+dxdt + - [ [ {v{x,t) - v^{x))~dxdt. 

'T Jo JSa ^0 JSa 

Prom (|6.3p . (16.4(1 and the passage to the limit as r ^ in (|6.52l) . we deduce that, 

- f A(t)dt ^0 as r ^ 0. 
^ Jo 



(6.52) 



(6.53) 



Thus, since A is a.e. non-increasing on (0,r), and A{t) > for a.e. t G (0, min(T, ^)), one 
then has 

A{t) = for a.e. t G (^0, min (^T, . 
Since a is arbitrary, we deduce that, 

u{x,t) < v{x,t) a.e. in Qt- 

□ 
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Remark 6.2 In llJj^ , the entropy process solution ^{x,t,a) was proved to be independent of a 
for a divergence-free function g G C^(Qr). However, for the case of a general non divergence- 
free function g £ W^^'°°(Qt); same result can be shown by adapting the same proof as in 
Theorem 3] taking into account the slight modifications that could be deduced from the proof of 
Theorem (2.16\) . More precisely, the treatment of the two terms and S'tep 2. 
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